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I A sample of edge-on spiral galaxies aimed at a thorough study of the main structural 

and photometric parameters of edge-on galaxies both of early and late types is pre- 
sented. The data were taken from the Two Micron All Sky Survey (2MASS) in the J, 
H and Kg filters. The sources were selected according to their angular size mainly on 
the basis of the 2MASS-selected Flat Galaxy Catalog (2MFGC). The sample consists 
of 175 galaxies in the if^-filter, 169 galaxies in the iJ-filter and 165 galaxies in the 
J-filter. We present bulge and disc decompositions of each galaxy image. All galaxies 
(~| I have been modelled with a Sersic bulge and exponential disc with the BUD DA v2.1 

Q^. package. Bulge and disc sizes, profile shapes, surface brightnesses are provided. 

Our sample is the biggest up-to-date sample of edge-on galaxies with derived 
structural parameters for discs and bulges. In this paper we present the general re- 
^ suits of the study of this sample. We determine several scaling relations for bulges 

' and discs which indicate a tight link between their formation and evolution. We show 

that galaxies with bulges fitted by the Sersic index n < 2 have quite different dis- 
tributions of their structural parameters than galaxies with n > 2 bulges. At a first 
approximation the Sersic index threshold n ~ 2 can be used to identify pscudobulgcs 
and classical bulges. Thus, the difference in parameter distributions and scaling re- 
lations for these subsamples suggests that two or more processes are responsible for 
I disk galaxy formation. The main conclusions of our general statistical analysis of the 

sample are: 

(1) The distribution of the apparent bulge axis ratio qh for the subsamplc with 
n < 2 can be attributed to triaxial, nearly prolate bulges that are seen from different 
projections, while n > 2 bulges seem to be oblate spheroids with moderate flattening. 

I Triaxiality of late-type bulges may be due to the presence of a bar that thickened in 

the vertical direction during its secular evolution. 

(2) For the sample galaxies, the effective radius of the bulge re,b, the disc scale- 
length h and the disc scaleheight zq are well correlated. However, there is a clear trend 
for the ratio r^.h/h to increase with n. As n is an indicator of the Hubble type, such a 

' trend unambiguously rules out the widely discussed hypothesis of a scale-free Hubble 

sequence. The found correlation between zq and r^b is new and was not described 
earlier. 

(3) There is a hint that the fundamental planes of discs, which links only disc 
parameters and the maximum rotational velocity of gas. are different for galaxies with 
different bulges. This may indicate a real difference of discs in galaxies with low and 
high concentration bulges. 

(4) The most surprising result arises from the investigation of the Photometric 
Plane of sample bulges. It turns that the plane is not flat and has a prominent curvature 
towards small values of n. For bulges this fact was not noticed earlier. 

(5) The clear relation between the flattening of stellar discs h/ zq and the relative 
mass of a spherical component, including a dark halo, is conflrmed not for bulgeless 
galaxies but for galaxies with massive bulges. 

Many of our results are in good agreement with the results of other authors, 
several ones are new. Thus, our sample is very useful for further detailed studying and 
modelling of the edge-on spiral galaxies. 

Key words: galaxies: fundamental parameters — galaxies: general — galaxies: in- 
frared — galaxies: spiral edge-on galaxies — galaxies: photometry. 
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1 INTRODUCTION 

The photometric and kinematic study of spiral galaxies pro- 
vides a method to tackle the problems concerning their for- 
mation and evolution. A lot of observational data has been 
obtained in recent decades, especially with the advent of 
CCD arrays. Over the last years, significant progress has 
been made in both observational and theoretical studies 
aimed at understanding the structure and evolution of galax- 
ies. Modern sky surveys (2MASS, SDSS, etc.) provide in- 
formation on the characteristics and spatial distribution of 
millions of galaxies. The large available amount of informa- 
tion has altered the face of modern extragalactic astronomy. 
Many important tasks, including detailed photometric stud- 
ies of the structure of nearby galaxies, can be now solved 
without additional observations by telescopes. Here we fo- 
cus on the surface photometry of spiral galaxies on the basis 
of the 2MASS near-infrared survey. 

The parametrization of galaxies is a way to describe 
the multitude of normal galaxies without any peculiarities 
within their structure. Various fitting functions are used 
to derive structural parameters of bulges and discs. Tra- 
ditional fitting functions for stellar discs include an expo- 
nential function (de Vaucouleurs 1959, Freeman 1970) and 
an inner-truncated exponential function (Kormendy 1977). 
Bulges are well described by the so-called Sersic profile that 
includes the de Vaucouleurs profile as a special case (Sersic 
1968). The de Vaucouleurs law, or r^^^-law, was obtained 
empirically for elliptical galaxies (de Vaucouleurs 1948) and 
later for bulges of early- type galaxies (de Vaucouleurs 1959). 
Later it was shown that the Sersic law is preferable to fit 
spheroidal components of galaxies (e.g. Andredakis et al. 
1995). The ^-distribution of surface brightness of edge-on 
discs is exponential or the so called "sech^-distribution" (van 
der Kruit & Searle 1981). (All the mentioned functions are 
described in Section 3 of this paper.) 

In order to derive structural parameters of bulges and 
discs, we must perform photometric bulge/disc decompo- 
sition of galaxy images. There are numerous methods of 
such decomposition. The one-dimensional analysis of the 
surface-brightness distribution is the most commonly used 
method (e.g. Boroson 1981; Kent 1985; Baggett, Baggett 
& Anderson 1998; Bizyaev & Mitronova 2002, hereafter 
BM02). But as was shown by Byun & Freeman (1995) the 
one-dimensional decomposition procedure leads to strong 
systematic errors. The two-dimensional approach provides 
much better estimates of the structural parameters for each 
component of a galaxy. This method is powerful when we 
need to separate the non-axisymmetric structures (bars, 
rings, spiral arms) from axisymmetric components (disc and 
bulge). 

In recent decades a number of studies concerning struc- 
tural properties of galaxies has been performed. The results 
of decompositions of near face-on disc dominated galax- 
ies as well as of early-type galaxies were presented, for in- 
stance, by de Jong (1996). It was concluded that there is 
a correlation between structural parameters of discs and 
bulges. MoUenhoff & Heidt (2001) found similar correlations 
and speculated about the formation and evolution of spi- 
ral galaxies. Graham (2001) performed a detailed analysis 
of low-inclination spiral galaxies and described many signifi- 
cant properties of their bulges and discs. Other studies dealt 



with the correlation between the scalelengths of bulges and 
discs (e.g. Courteau et al. 1996; de Jong 1996), the presence 
of additional morphological components such as bars, lenses, 
rings, inner bars and inner discs (e.g. de Souza, Gadotti 
& dos Anjos 2004, SGA04 hereafter; Castro-Rodn'guez & 
Garzon 2003; Laurikainen, Salo & Buta 2005; Gadotti 2009) 
and the revision of morphological classification of galaxies 
(SGA04). The fundamental plane (FP) of discs as well as 
the photometric plane of bulges were also widely discussed 
(e.g. Khosroshahi et al. 2000a,b; Mollenhoff & Heidt 2001; 
Moriondo, Giovanelli & Haynes 1999). 

The main aim of our project is to perform decompo- 
sitions of a large number of edge-on galaxies both of early 
and late types into bulges and discs. We chose this type of 
galaxies because they provide a unique possibility to obtain 
information about the vertical structure of galactic discs. In 
the edge-on view non-axisymmetric features such as bars, 
rings or spiral arms are not seen. On the other hand, the 
vertical distribution of light in edge-on galaxies allows us 
to study a thin disc, a thick disc and a stellar halo (see 
e.g. Yoachim & Dalcanton 2005; Seth, Dalcanton & de Jong 
2005). 

Unfortunately, there are some difficulties in studying 
edge-on galaxies. In the optical band we need to take into 
account the dust extinction in the plane of a disc. The distri- 
bution of dust is often unknown and it is hard to obtain true 
values of disc and bulge parameters. The advantage of NIR 
observations is that the general extinction is drastically re- 
duced, although we can not exclude it completely. Neverthe- 
less it is better to perform the decomposition of galaxies by 
using their images in red and infrared bands. Using three IR 
bands for galaxy decomposition and subsequent statistical 
analysis of galaxy properties in all three bands encourages 
us in our suggestion that extinction problems are not very 
severe (at least for general statistical correlations which are 
the subject of the present work). However, edge-on galax- 
ies are rare in occurrence in contrast to face-on galaxies or 
galaxies with intermediate inclination angles. Therefore, in 
order to create a large sample of edge-on galaxies it is nec- 
essary to use special catalogues which include a lot of edge- 
on spiral galaxies such as the Revised Flat Galaxy Catalog 
(Karachentsev et al. 1999). 

A detailed study of the disks of edge-on galaxies is de- 
scribed in many papers (e.g., van der Kruit & Searle 1981; 
Reshetnikov & Combes 1997; de Grijs 1998). In a more re- 
cent study, Dalcanton & Bernstein (2000, 2002) examined a 
sample of extremely late- type, edge-on spirals with no ap- 
parent bulges and concluded that thick disks are common 
around galaxy disks of all masses. BM02 also analyzed a 
sample of late-type edge-on galaxies in the J, H and Ks- 
bands to compare the mean ratios h/zo in these bands. They 
have noted a strong correlation between the central surface 
brightness of the disk and the h/zo ratio: the thinner the 
galaxy, the lower the central surface brightness reduced to 
the face-on inclination. Kregel, van der Kruit & de Grijs 
(2002), Kregel (2003), Kregel, van der Kruit & Freeman 
(2005) and Zasov et al. (2002) presented a complex study 
(dynamical properties of the stellar discs, three-dimensional 
disc structure, stellar kinematics, rotation curves) of several 
dozens of edge-on spirals. Unfortunately, there are no ex- 
tended studies of structural parameters of edge-on galaxies 
both of early and late types which join analysis of their 
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bulges and discs. Our study is intended to fill this gap. 
We analyse the statistical properties of the sample edge-on 
galaxies and present correlations between the global struc- 
tural parameters of bulges and discs, including those that 
describe their vertical structure, in the near-infrared bands. 

This paper is organized as follows. In Section 2 we 
present the detailed description of our sample of spiral edge- 
on galaxies which includes the selection criteria and a brief 
description of the data sources. We also describe the com- 
pleteness of the sample and its main properties. In Section 3 
we use the two-dimensional least-squares algorithm to de- 
compose the galaxy images into a bulge and a disc. We recall 
some aspects of the code BUDDA (SGA04) which we used 
as the decomposition tool. In Section 4, we present the re- 
sults of image decomposition in the J, H and Js's-bands. We 
compare our results with previous results from the litera- 
ture and analyse the reliability of the results of the BUDDA 
decompositions. In Section 5 we describe general properties 
of edge-on spiral galaxies and discuss some correlations be- 
tween galactic subsystems. Some of them are well-known 
but there are also some new results. In Section 6 we sum- 
marize our main conclusions. In the forthcoming papers we 
are going to discuss properties of edge-on galaxies in more 
detail. 



2 THE SAMPLE 

We intended to find edge-on galaxies with a wide range of 
bulges: from bulge-dominated spirals to galaxies with small 
bulges. We used the 2MASS-selected Flat Galaxy Catalog 
(2MFGC) (Mitronova et al. 2003) as a source of objects and 
the 2MASS survey (Skrutskie et al. 2006) as a source of 
image data. We also used the Revised Flat Galaxies Cata- 
log (RFGC) by Karachentsev et al. (1999) to add late-type 
galaxies to our sample. 

The Two Micron All Sky Survey is one of the best 
known digital surveys in a wavelength range close to the 
optical one. This survey covers the whole sky in the filters J 
(1.25 ^lm), H (1.65 ^m), and Ks (2.16 fim). Unfortunately, 
the 2MASS survey has a weak sensitivity to late-type galax- 
ies, especially of low surface brightness (Jarrett 2000). This 
is due to the high brightness of the night sky in the near-IR 
range and short exposures. For this reason the periphery of 
the discs of spiral galaxies is generally unseen beyond the 
isophotes fainter than = 20 mag arcsec"^. We did not 
study truncation radii for the 2MASS-galaxies and did not 
take them into account performing the decomposition. 

The sample was selected mainly from the 2MFGC. This 
catalogue contains 18020 disc-like galaxies covering all the 
celestial sphere. The objects were taken from the Extended 
Source Catalog of the 2MASS (XSC 2MASS) according to 
their 2MASS axial ratio a/6 ^ 3. The 2MFGC gives value of 
axial ratio h/a, fiducial elliptical Kron magnitude, position 
angle, the J-band concentration index (/C) and other useful 
information. Concentration index is the ratio of the radii of 
circles that contain 3/4 and 1/4 of the measured flux in the 
same band, respectively. IC together with b/a are known to 
serve as dividing lines between early- and late-type galaxies 
(see e.g. Strateva et al. 2001; Nakamura et al. 2003). We 
decided to select galaxies both of early and late types in 
almost equal quantities. Our criteria are 



- axial ratio for the combined J+H + K image, sba > 0.2; 

- if-band fiducial elliptical Kron radius, Rk ^ 30"; 

- Hubble type ranging from SO to late types; 

- non-interacting galaxies; 

- concentration index IC > 2.0. 

The latter item means that we exclude galaxies with 
strongly peculiar surface brightness distributions since for 
pure exponential disks IC = 2.8 and for r^^* ellipticals IC — 
7 (Fraser 1972). 

From the set of 3167 objects which contains galaxies 
with pronounced bulges (Sa, Sb types) and galaxies with 
small bulges (Sc, Scd types) we selected a sufficient num- 
ber of regular spiral galaxies to make statistical analysis of 
edge-on galaxies properties. We did not add to the sample 
galaxies that are clearly lopsided or warped (with some ex- 
ceptions below). In addition, on the basis of careful visual in- 
spection we rejected non-edge-on and peculiar objects. Also, 
we added 13 bright galaxies that were absent in the 2MFGC 
and several edge-on galaxies which were not selected accord- 
ing with foregoing criteria but were included in the sample 
of RFGC galaxies by BM02. We chose them to compare the 
results of our decomposition with those obtained by BM02. 

The final sample consists of 67 early-type galaxies (SO- 
Sab), 61 late-type galaxies (Sbc-Sd) and 47 galaxies of in- 
termediate type Sb. The absolute numbers of galaxies of 
various morphological types (taken from the LEDA) and 
their percentages are listed in the Table [1] The full list of 
the sample galaxies studied here is presented in the Ap- 
pendix, where we provide the basic information about them. 
According to LEDA 23 galaxies have bars, 6 galaxies have 
rings, and 3 objects are interacting galaxies. We did not take 
into account additional components like bars performing the 
decomposition, but checking all correlations between struc- 
tural parameters, we did not find any strong deviation of 
the above-mentioned objects from derived trends. 

We also did not reject 3 interacting galaxies and several 
other galaxies with some unusual structural features (double 
exponential discs, rings, etc., which will be described in our 
resulting tables) to have about 10 unusual spiral edge-on 
galaxies for further work. 

As can be seen in the Table[T]the sample contains a sub- 
stantial percentage of lenticular galaxies. It is necessary to 
notice that Hubble classification for edge-on galaxies is very 
subjective because the spiral arm structure is not seen in 
this case. The only way to classify edge-on galaxies accord- 
ing to morphological types is the estimation of the ratio of 
the total luminosities of a bulge and a disc (hereafter B/D). 
We assume that a large fraction of classified SO galaxies are 
Sa galaxies in fact. In Fig. [T] the images of typical sample 
galaxies are shown. According to the LEDA data, our sample 
galaxies are, on average, true highly-inclined galaxies with 
the mean inclination i = 87. °8 ± 4.°8. Of course, this esti- 
mate can be too optimistic but a moderate deviation from 
edge-on orientation does not significantly change the slope 
of the vertical surface brightness distribution (Barteldrees 
& Dettmar 1994). 

The distances of the galaxies are taken from 
NASA/IPAC Extragalactic Database (NED) with the Hub- 
ble constant Ho = 73 km s~^/Mpc, nmatter = 0.27, and 
f^vacuum ~ 0.73). The luminosity distances are in the range 
from 4.3 Mpc (NGC 891) to 171 Mpc (ESO 251-G028) with 
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Table 1. Morphological types of the 
sample galaxies and their fractions 



Type 


Number 


Percenta 


SO, SO-a 


41 


23.4 


Sa 


11 


6.3 


Sab 


15 


8.6 


Sb 


47 


26.9 


Sbc 


21 


12.0 


Sc 


36 


20.5 


Scd 


4 


2.3 


Total 


175 


100 




Figure 1. Examples of the sample galaxies in the Ks-ha.nd. The 
upper image is a typical galaxy of SO-a type: ESO 311-G012. The 
middle image is NGC 4013 of Sb type. The bottom image is an 
example of Sc galaxy: NGC 5907. 



a sample median value of 40 Mpc. Three galaxies (ESO 013- 
G024, ESO 555-G023, ESO 555-G032) have no measured or 
estimated distances. We do not take them into account in 
statistical analysis where the distances are required. 

To investigate some kinematical and dynamical prop- 
erties of the sample galaxies we added information about 
maximum rotation velocity for each galaxy if it was pre- 
sented in the LEDA data base (see the section 5.3 for more 
details). 

The sample consists of 175 galaxies in the Jfs-filter, 169 
galaxies in the H-filter and 165 galaxies in the J-filter. (Sev- 
eral galaxies have prominent dust lanes within their discs 
and even in the J and H passbands absorption of the dust 
is significant.) 

In order to evaluate the completeness of the sample, we 
use the method proposed by Schmidt (1968) (the so-called 
V/Vmax completness test). Following Thuan & Seitzer (1979) 
we used this method for our angular-size-limited sample. 
The full sample is incomplete but the subsample of galaxies 
with angular radius r ^ 60" is complete (1^/Vmax = 0.49 ± 
0.03). This subsample consists of 92 galaxies (47 early-type 
galaxies and 45 late- type galaxies). 

The images used in this work were taken from the 
database of the All-Sky Release Survey Atlas. The Atlas' 
images are FITS standard images in 512 x 1024 pixel format 
(with 1. 00 pixeP^). Description of the images such as the 
pixel gain, readout noise, seeing FWHM (~ 2. "5), sky value 
and also zero points are available in the 2MASS webpages. 
This information has been used in making decompositions 
of galaxy images. Each image was thoroughly examined in 
order to mask background stars and galaxies and other con- 
taminating sources. All the images were rotated to align the 
galaxy major axis parallel to the horizontal image borders. 



These procedures were carried out in three filters so all im- 
ages of the same galaxy in the J, H and Ks bands have the 
same extension and position angles. The smallest image of 
our sample is 48 x 18 pixels for ESO 555-G032 and the largest 
image is 869x165 pixels for NGC 4565. The median value of 
the size of galaxy images is 119x27 pixels, that is enough to 
determine the main photometric parameters of galaxy discs 
and bulges. The forthcoming analysis was performed with 
the MIDAS (developed by the European Southern Observa- 
tory) and the BUDDA packages. 



3 TWO-DIMENSIONAL BULGE/DISC 
DECOMPOSITION 

BUDDA (Bulge/Disc Decomposition Analysis) is a code de- 
vised to perform two-dimensional decomposition of galac- 
tic images. This code was produced by R.E. de Souza, 
D.A. Gadotti and S. dos Anjos (for more information about 
BUDDA see SGA04). Since the year of 2004 BUDDA has 
been publicly available to the astronomical communitjQ. We 
use BUDDA v2.1 to produce the fits, allowing the inclusion 
of bulges, discs, bars and AGN in the models. 

The input data for the code is an image of a galaxy, 
consisting of two major components: a bulge and a disc. 
The disc is represented by an exponential distribution of the 
luminosity density J(r, z) with the central luminosity density 
/(O, 0), the scalelength h and the 'isothermal' scaleheight zq: 

lir, z) = J(0, 0) e-"/'^ sech\z/zo) , (1) 

where (r, z) are the cylindrical coordinates. The disc is as- 
sumed to be axisymmetric and transparent. 

The surface brightness distribution for the face-on disk 
(i = 0°) can be expressed as follows: 

/d(o = e"'■''^ (2) 

where Jo = 2 zo 1(0,0). The same expression in magnitudes 
per arcsec^ is: 

Atd(r) = So,d-f 1.086 r-//i, (3) 

where So.d is the central surface brightness of the face-on 
disc. For edge-on galaxies (i = 90°) the following expression 
is valid: 

7(r,z) =1(0,0)^7^1 (^) sech\z/zo), (4) 

where K\ is the modified Bessel function of the first order. 

The bulge surface brightness profile is described by the 
Sersic law (Sersic 1968): 

/b(r) = /o,be-"^"'^"''''-^>''"', (5) 

where r^^h is the effective radius of the bulge, i.e., the radius 
of a circle that contains 50% of the total galaxy luminosity, 
7o,b is the bulge central surface brightness, n is the Sersic 
index, defining the shape of the profile, and the parameter 
Vn ensures that re,b is the half-light radius. In magnitudes 
per arcsec^ the expression looks as follows: 

Mb(r) = Mo,b + 3^ — j , (6) 
^ see |http://www.mpa-garching.mpg.de/'^dimitri/budda.html| 
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where /^o,b is the bulge effective surface brightness ex- 
pressed in mag per arcsec'^, i.e., the surface brightness at 
re,b. The BUDDA code uses a numerical approximation of 
u„ ~ (0.868ri - 0.142) In 10. The Sersic index of n = 4 repre- 
sents the de Vaucouleurs profile which was very popular to 
describe the surface brightness distribution of bright ellip- 
tical galaxies and of bulges in early-type spirals. The Sersic 
index of n = 1 represents the exponential profile of bulges 
in late-type spirals, of galactic discs and of dwarf elliptical 
galaxies (see e.g. de Vaucouleurs 1948, 1959; Freeman 1970; 
Graham 2001 and references therein). 

It is also useful to determine such parameters as ellip- 
ticities (e = 1 — b/a) of a bulge eb and a disc ed, as well as 
their position angles P.A.b and P.A.^- To account for the 
effect of seeing there is a Moffat smearing (Moffat 1969) of 
the brightness profiles in the code controlled by a parameter 
associated with the seeing radius as. There is also included a 
correction term of the sky level Aaky The isophotes of bulges 
and discs can be described by generalized ellipses with the 
ellipse index parameter 7 controlling their shape: 



Table 2. Parameters of the two-component model for 
ESO 240-GOll (Ks band) 



= 1, 



(7) 



where x and y are the pixel coordinates of the ellipse points, 
b and a are the sizes of its semi-major and semi-minor axes, 
respectively. When 7=1 the ellipse is simple, while for 
7 > 1 the ellipse is boxy and for 7 < 1 the ellipse is discy. 

For edge-on galaxies BUDDA derives the disc scale- 
height zo instead of the disc ellipticity e^- Thus, a total 
number of the main parameters for the bulge/disc model 
is 11: the center coordinates xo and yo, the central surface 
brightnesses of a bulge /o,b (or the effective surface bright- 
ness of bulge 7e,b as used in BUDDA) and a disc /o,d (both 
in ADU), rc,b, h, P.A.h, P-A.d, the disc scaleheight zq, the 
bulge ellipticity eb and the Sersic index of a bulge n. 

Total luminosity of a disc can be expressed as: 



27r Jo.d h 



while the total luminosity of a bulge is 

= r(27i) Jo,b re,b (1 - Eb) 



(8) 



(9) 



where F means the gamma function. 

The code BUDDA allows the user to add some addi- 
tional components in a galaxy model (bars, AGNs and dou- 
ble exponential discs), but we made decompositions of all 
our galaxies into a bulge and a disc only. There are some 
difficulties to determine the presence of a bar or any other 
substructures in edge-on perspective. The observed surface 
brightness distribution is a result of luminosity integration 
along the line of sight and such integration can smooth and 
hide such subcomponents. It is also hard to make the first 
approximation of the values of parameters which character- 
ize these additional features. 

Fig. [2]illustrates typical results of our decomposition for 
the Sb-type galaxy ESO 240-GOll. The size of the image is 
228x40 pixels. The sky background. Galactic stars, bright 
HII regions and nearby companions were removed. The de- 
rived parameters are given in Table (2] As seen in Fig. [2] the 
fit is fairly good. 



D (Mpc) 35.7 



Me.b 



fe,b (kpc) 

^J.O,d 16.22 ±0.16 n 

h (kpc) 5.1 ± 0.9 (?b 

zo (kpc) 0.80 ±0.13 B/D 

h/zo 6.4 Drot (km/s) 



17.09 ± 0.09 
1.34 ±0.09 
2.4 ±0.3 
0.69 ±0.04 
0.38 

267.5 ± 3.3 



4 DECOMPOSITION AND RESULTS 

To start the galaxy decomposition, we must determine the 
first-order approximations for the input parameters (see pre- 
vious section). To obtain the first-order approximation for 
our 2D fits, we have used parameters derived from ID fits 
of the major and minor axes profiles. 

In almost all cases the ellipse index 7 of a bulge was 
taken to be equal 1.0 ± 0.01 (for simplicity we consider 
bulges with elliptical isophotes). A galaxy with complex 
morphology of a bulge requires more detailed investigation 
with higher resolution images than those that are given by 
2MASS. The result of decompositions into a bulge and a disc 
is the list of all structural parameters of each object of the 
sample in 3 photometric bands. The final values of ^c,h and 
Po,d were corrected for the Galactic extinction according to 
Sciilegel, Finkbeiner & Davis (1998). 

To check the reliability of our 2D decompositions, we 
subtracted the full apparent magnitudes of our models from 
the apparent magnitudes presented in the 2MASS catalogue. 
The mean values and standard deviations of these subtrac- 
tions for the sample galaxies are: A J = 0."'011 ± O.^OS?, 
AH = 0.^021 ± 0.^049 , AKs = 0.^014 ± 0.™053 . 

Our sample has 30 galaxies in common with the sample 
of BM02. The method used by these authors is based on 
several photometric cuts parallel to the minor axis of each 
galaxy and parallel to its major axis. Their sample includes 
flat galaxies mainly of late types in the Ks filter. Fig.Oshows 
a comparison of radial scalelengths (the upper frame), verti- 
cal scale heights (the middle frame) and deprojected central 
surface brightness of discs (the bottom frame). Open squares 
are related to the objects which have some peculiarities or 
features in structure (NGC 5965, UGC 6012, UGC 12533 are 
with possible bars, IC 4202, UGC 4517, ESO 121-G006 are 
with strong dust lanes, and UGC 1817 is a bulgeless galaxy 
in a group). These objects do not lie on the diagonal. The 
mean values and standard deviations of difference between 
parameters obtained by us and by BM02 are: 

(/i-Zibm) = 2."46±2."11, 
(zo -^o,bm) =0."41±0."29, 

(50 -So.bm) = 0.^28 ±0.^19. 

As one can see (Fig. |3} , our results and those of BM02 are 
in good agreement. 



De Grijs (1998) presented detailed surface photometry 
for a sample of edge-on galaxies in the B, I and K bands. 
Fig- 51 shows a comparison of radial scalelengths for 14 joint 
galaxies (Tv-band). Some systematics are evident which can 
be attributed to different decomposition procedures. 
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r (arcsec) z (arcsec) 

Figure 2. An example of the edge-on spiral galaxy ESO 240-GOll in the i^s-band. The images on the top show the galaxy (top left), 
the model (top middle) and the residual image (top right). The residua were obtained after subtraction of the model from the image. 
The bottom panels show the surface brightness profiles of the galaxy, the model and the residual images arc indicated. The profiles are 
drawn along major- (left) and minor axes (right) at r = 0.5h, l.Oh and l.5h. 



Table [3] compares decomposition results for 5 galaxies 
in the present work and in the paper by Bedregal, Aragon- 
Salamanca & Merrifield (2006). As one can see, the results 
are in agreement that is typical for such a comparisoijfj. It 
proves the reliability of the obtained results and makes them 
suitable in our further statistical analysis. 



5 MAIN RESULTS 

Statistical correlations between structural parameters give 
us a clue for understanding the formation and evolution of 
spiral galaxies. Here we present only some of the correla- 
tions of the main parameters of the galaxies. The detailed 
analysis of the properties of bulges and discs will be done in 
a forthcoming article. 

5.1 Distributions of structural parameters 

Here we briefly describe the distributions of structural pa- 
rameters of stellar components. In general, we do not take 
into account any corrections for the incompleteness and se- 
lection effects of the sample but in special cases we give the 
best fitting parameters obtained for the complete subsam- 
ple. 

5.1.1 Discs 

The distributions of the apparent central surface brightness 
of discs for edge-on view /io,d corrected for Galactic extinc- 
tion in the JHKs filters are shown in Fig. [5^. The values of 



this parameter span more than 3™ arcsec"'^ in each filter. 
The distributions shift from the filter J to the filter Ks with 
increasing of mean central surface brightness. It is well seen 
that LSB galaxies were not included in the sample. Thus, 
one has to keep in mind that in fact the number of such 
galaxies could be high at the faint end of the distribution of 
Ho.d in all three passbands. 

Fig. [5}3 shows a large range in scalelengths h. Most of 
the sample galaxies have h — 1 — 6 kpc but the number 
of extended galaxies is small. The distributions of the disc 
scaleheight zo are shown in Fig. [5j;. There is a fairly narrow 
distribution, that lies in the range zq = 0.4 — 1.2 kpc, and 
a tail of "thick" galaxies with zq > 1.2 kpc. The shapes of 
distributions of h and zo in each filter are quite similar and 
the median values of h and zo have no systematic differences 
in the infrared bands under consideration (Table |4]). 

The statistical study of the ratio of the radial to the ver- 
tical scalelengths for galactic discs, h/zo, gives constraints 
on kinematical and dynamical models of galaxies. The dis- 
tribution of this ratio h/zo is shown in Fig. [5]i. It is rather 
flat. The mean ratio (h/zo) is about 3.9 for the H and Ka- 
bands (Table[B)). Kregel et al. (2002) found the mean ratio of 
{h/zo) = 3.7±1.1 (in the /-band) which is close to our result. 
For the sample of 153 galaxies composed by BM02 the mean 
ratio of h/zo is about 4.8. Although they studied galaxies 
in the ifs-band, their sample consists of late-type galaxies 
and the mean value of the ratio h/zo for these galaxies is 
larger than the mean ratio found in this paper. Apparently 
the mean h/zo ratio is higher in the blue filters than in the 
near-infrared bands (see for comparison e.g. van der Kruit 
& Searle 1982a and de Grijs & van der Kruit 1996). 



2 For NGC 5308 and NGC 1184 Bedregal et al. (2006) give the 
formal errors for the bulge effective radii that are very large and 
asymmetric. Unfortunately, they did not provide any explanation 
of these facts. For both galaxies, the bulge effective radii from our 
study and those by Bedregal et al. (2006) are quite different. 



5.1.2 Bulges 

The distributions of the apparent effective surface brightness 
of bulges are shown in Fig. [5^. The distributions in each filter 
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Table 3. Comparison between structural parameters obtained in this work and by Bedregal ct al. 
(2006) in the Ka-hand (parameters given in arcsec were rewritten in lipc according to distances 
adopted by Bedregal ct al.). 
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Figure 3. Comparison of the radial scalelengths (the upper 
frame), the vertical scaleheights (the middle frame) and the de- 
projected central surface brightnesses (the bottom frame) be- 
tween data of BM02 and the data presented in this paper. Open 
squares are related to the galaxies with some peculiarities and 
features within their structure as noted in Sect. 4. 



Figure 4. Comparison between radial scalelengths obtained by 
de Grijs (1998) and in this work. 



have a narrow peak in contrast to the distributions of the 
central surface brightness of discs. 

Fig. [5f demonstrates the fairly broad distributions of the 
effective radii of the bulges of galaxies. The distributions of 
the Sersic indices of galaxies are shown in Fig. [S^. The dis- 
tributions demonstrate a weak bimodality which may reflect 
the existence of two families of bulges: bulges with n > 2 
and bulges with n < 2. This bimodality is known for some 
galaxy samples (e.g. Fisher & Drory 2008) and probably is 
real for our sample. As was reported earlier (see e.g. Fisher 
& Drory 2007, 2008) such a bimodality correlates with mor- 
phological type of the bulge. Classical bulges have n > 2 
and so-called pseudobulges have n < 2, but this Sersic in- 
dex threshold can be considered only as an approximation to 
identify pseudobulges. Further we are showing the difference 
in the photometric planes for these two types of bulges. 

In Fig. [SJi we demonstrate the distributions of model 
bulge axis ratio gb = for edge-on view. In the case of 
edge-on galaxies this parameter describes directly the bulge 
flatness and the intrinsic 3D structure of bulges if they are 
assumed to be oblate spheroids. The Table [6] gives the me- 
dian value of this parameter ~ 0.63, independently of the 
band. This is in good agreement with Moriondo, Giova- 
nardi & Hunt (1998) and Noordermeer & van der Hulst 
(2007). Moriondo et al. (1998) derived for their samples of 
14 moderately inclined galaxies the median intrinsic axis ra- 
tio gb = 0.64 while Noordermeer & van der Hulst (2007) 
gave the average value of (gb) = 0.55 ±0.12 for their sample 
of 21 early-type disc galaxies decomposed into contributions 
from a spheroidal bulge with a Sersic proflle and a flat disk. 
Thus, bulges are definitely nonspherical and flattened. 

In recent years a great deal of new data has made 
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researchers revise the assumption that bulges are oblate 
spheroids. However, conclusions about intrinsic 3D shapes of 
bulges and their possible triaxiality are still somewhat con- 
troversial, despite the fact that such conclusions are thought 
to be crucial for testing different scenarios of galaxy for- 
mation. For galaxies at intermediate inclination angles mis- 
alignment between the bulge and the disk major axes in- 
dicates that bulges are probably triaxial (Bertola, Vietri 
& Zeiliger 1991). From the observed distribution function 
of apparent ellipticities and misalignment angles Bertola et 
al. (1991) determined the probability distribution function 
(PDF) of intrinsic axial ratios for 32 bulges of their sample. 
The peak of the PDF falls on the intrinsic equatorial axis 
raticO B/A = 0.95 and bulge flattening C/A = 0.65. About 
one half of all bulges seems to be close to oblate, with the re- 
mainder being triaxial. Fathi and Peletier (2003) measured 
the bulge deprojected ellipticity in the equatorial plane that 
was obtained from the ellipse fitting the galaxy isophotes 
within a bulge region. One can estimate the median of B/A 
from their distribution function as 0.7 — 0.8 with smaller val- 
ues attributed to more elongated bulges in late-type galax- 
ies. Mendez-Abreu et al. (2008) also focused their atten- 
tion on the intrinsic equatorial ellipticity of bulges. They re- 
constructed the PDF for the sample of 148 unbarred SO-Sb 
galaxies and obtained the value of {B/A) = 0.85. This is con- 
sistent with findings by Bertola et al. (1991) and Fathi and 
Peletier (2003). However, in contrast to Fathi and Peletier 
(2003), Mendez-Abreu et al. (2008) did not find any signif- 
icant differences in the shape of bulges between samples of 
early- and late- type galaxies. 

Our sample of edge-on galaxies allows us to distinguish 
clearly the difference in the bulge equatorial ellipticities for 
early- and late-type galaxies. As a Sersic index is correlated 
with the Hubble type (Andredakis, Peletier & Balcells 1995; 
Graham 2001; MoUenhof 2004) (as well as with the bulge 
morphology, dividing bulges into classical ones and pseudob- 
ulges; Fisher & Drory 2008) we divided our sample bulges 
between those with Sersic index n > 2 (early-type galaxies 
or classical bulges) and those with n < 2 (late-type galaxies 
or pseudobulges) using Ks filter. As one can see in Fig. |6] for 
bulges with n > 2 the distribution of qh has a rather narrow 
peak at qh ~ 0.65. This may reflect the fact that bulges in 
early-type spirals are nearly oblate spheroids with moderate 
flattening. The distribution of qh for bulges with n < 2 is 
very wide, spreading from flat bulges up to nearly spheri- 
cal ones. Such a distribution may be attributed to definitely 
triaxial, near prolate bulges that are seen from different pro- 
jections — along the major axis and perpendicular to it. But 
we can not exclude that triaxial shape of bulges in late-type 
galaxies may hide the presence of bars that thickened in the 
vertical direction during secular evolution. We are planning 
to investigate this question in more detail in a forthcoming 
paper. 

Fig. [7] demonstrates the distributions of the ratio of 
bulge and disc luminosities B/ D = Lb/Ld- Our sample con- 
tains many early-type galaxies with B/D > 0.5 and late- 

^ We use capital letters A, _B, C to define the intrinsic semi-axes 
of a triaxial system and small letters a, b for the apparent semi- 
axes of an image. In axisymmetric case and an edge-on view the 
ratio bja coincides with the ratio C/A. 



60 




0.2 0.4 0.6 0.8 1 



Figure 6. Distribution of the sample galaxies over the model 
bulge axis ratio qh in the i^s-band. The thin line indicates the 
distribution for the whole sample, the solid line corresponds to 
the subsample of bulges with n > 2 and the dashed line shows 
the distribution for the subsample of bulges with n < 2. 



type galaxies with values of B/D < 0.2. For edge-on galaxies 
the ratio of B / D is a key parameter in addition to n for clas- 
sifying spiral galaxies on the Hubble sequence because the 
spiral arms are not seen in this case (de Jong 1996; Graham 
2001). Our sample also demonstrates a correlation between 
the ratios oi B/D and the Hubble morphological types T 
taken from LEDA (Fig. [SK) that is in agreement with pre- 
vious work (e.g., de Jong 1996; Graham 2001; MoUenhoff 
& Held 2001; Hunt, Pierini & Giovanardi 2004; MoUenhoff 
2004) . However for edge-on galaxies morphological types are 
very subjective. Comparisons of the Hubble types given by 
different classifiers show an rms uncertainty in the type in- 
dex of order 2 T- units (Lahav et al. 1995). Unfortunately 
there are errors in the bulge/disc decomposition. Therefore 
we use both approaches (type T from LEDA and obtained 
B/D) for describing morphological types of edge-on galax- 
ies. In some cases we will distinguish early- and late-type 
galaxies using the Sersic index n because there is a good 
correlation between B/D and n (Fig. [SJd). This correlation 
is not new and was presented by many authors (e.g., An- 
dredakis et al. 1995; Graham 2001; MoUenhoff & Held 2001; 
Hunt et al. 2004). 

Details of the statistical properties of the sample discs 
and bulges are presented in the Table lU the Tableland the 
Table |6l The values of n in Table [5] systematically increase 
going from J to Ks - this may imply that n is systematically 
underestimated due to extinction affecting the central peak 
of emission at the shorter near-IR wavelengths. Another pos- 
sible explanation is a real change of the surface brightness 
shape with wavelength. But, as one can see in Table (5] the 
quoted changes are within statistical errors. 

5.2 Scaling relations 

The relationships between bulge and disc parameters might 
give an insight in the ages of the bulge and disc formation 
and help to reconstruct the chronology of these events. An- 
dredakis & Sanders (1994) were the first to report the ex- 
istence of a correlation between disc and bulge scalelengths 
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Figure 5. (a) Distributions of the sample galaxies over (a) the central surface brightness of the disc, (b) the disc scalelength, (c) the 
disc scaleheight, (d) the h/zo ratio, (e) the effective surface brightness of the bulge, (f) the bulge effective radius, (g) the Sersic index 
of bulge, (h) the model bulge axis ratio gi, = b/a. In each plot the solid line corresponds to the Kg-hand, the dashed line shows the 
distributions corresponds to the //-band and the dotted line corresponds to the J-band. 
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Figure 7. Distribution of the sample galaxies over the ratio 
bulgo-to-disc luminosities B/D = L^^/L^. 



Table 4. The median values of ^iq zq and Af(j and their 

quartiles (75 and 25 per cent, respectively) 
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for their sample of 34 late-type spirals. Just after that, such 
a correlation was confirmed and used for speculating about 
the scale-free Hubble sequence since the relative size of bulge 
and disc seemed not to depend on the morphological type 
(de Jong 1996; Courteau et al. 1996). De Jong (1996) found 



Table 5. The median values of pL^ i,i '"c,bi ^ E'ld ^^b a-nd their 
quartiles (75 and 25 per cent, respectively) 
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Table 6. The median values of h/zp, (?b and 
B/D and their quartiles (75 and 25 per cent, 
respectively) 
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the equation for the least squares fitted line in the K pass- 
band: lg(rc,b) = 0.95 lg(/i) — 0.86 with a standard deviation 
of 0.17 and correlation coefficient r ~ 0.8. This suggests a 
linear relation between h and r^.h, i.e. the disc scalelength 
increases when the bulge effective radius increases. 

Furthermore it became clear that the ratio r^^h/h ap- 
peared to be strongly infiuenced by the bulge parametriza- 
tion (Graham & Prieto 1999; Graham 2001). The indepen- 
dence of morphological type and a strong correlation were 
seen only when galaxies were best modeled with a fixed 
n = 1 bulges for all type galaxies (de Jong 1996; Courteau et 
al. 1996; Moriondo et. al. 1998; Hunt et al. 2004). According 
to whether bulges were fitted with the n = 1 Sersic model or 
with the de Vaucoulers law, the correlation between scale- 
lengths was clear or disappeared completely (Andredakis & 
Sanders 1994; de Jong 1996; Hunt et al. 2004). 

Using Sersic bulges with a fi:ee shape parameter n as the 
best fitting models (MoUenhoff & Heidt 2001; Khosroshahi 
et al. 2000b; Graham 2001; MacArthur, Courteau & Hohz- 
man 2003) or different but fixed n for early- and late-type 
galaxies (Graham & Prieto 1999) confirms this correlation 
with a fairly large scatter and opens the discussion about the 
variation of the ratio re,b//i with Hubble types (Moriondo et 
al. 1998; Graham & Prieto 1999; Khosroshahi et al. 2000b; 
Graham 2001; Scodeggio et al. 2002; MacArthur et al. 2003; 
MoUenhoff 2004) and about the prematurity of the claims 
made by de Jong (1996) and Courteau et al. (1996) for a 
scale-free Hubble sequence. 

We also found that there exist fairly strong correlations 
between structural parameters of discs and bulges. The cor- 
relation between h and re,b is shown in Fig. |9^. In this figure 
and hereafter we mark the complete subsample galaxies with 
filled squares. The regression line is plotted only for points 
of the complete subsample. The equations of the linear re- 
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gression and the correlation coefficients are presented in the 
Table[7] These equations are consistent with a linear relation 
between h and re,b- 

We found a mean rc.h/h in the J-band of 0.39 ± 0.18, 
in the H-hanA of 0.42 ± 0.20 and in the ifs-band of 0.40 ± 
0.20, which are in agreement within errors with the values 
obtained by Mendez-Abreu et al. (2008) for their sample of 
148 unbarred SO-Sb galaxies (rc,h/h) = 0.36 ± 0.17 (in the 
J-band). 

There is some controversial evidence for systematic 
changes with morphological type in the ratio re,h/h (Scodeg- 
gio et al. 2002; MacArthur et al. 2003; Hunt et al. 2004; 
MoUenhoff 2004). Despite being difficult to detect, the de- 
pendence with Hubble types might be real in the sense that 
the sample of predominantly late-type disk galaxies gives 
{rs,h/h) = 0.22 ±0.09 (MacArthur et al. 2003; the _ff-band) 
and the sample of predominantly early-type galaxies leads 
to {rc,b//i) = 0.33 ± 0.17 (Khosroshahi et al. 2000b; the K- 
band) . 

For edge-on galaxies the reliable indicator of morpholog- 
ical type is B/D or n. Fig.^jp demonstrates the dependence 
of the ratio Ve^h/h on the Sersic shape parameter n. There 
is a clear trend for r^^h/h to increase with n. The size ratio 
increases from n ~ 0.5 to ti ~ 5 by a factor of ~ 3 — 4. 
Almost all bulges of the sample of MacArthur et al. (2003) 
have n < 1.2 and most of the bulges of the sample of Khos- 
roshahi et al. (2000b) have n > 2.2. Thus, the difference 
in the mean value of re^b/Zi for these two samples reflects 
the same tendency that one can see in Fig. The same 
increase was declared by Hunt et al. (2004). They found 
for n = 1 bulges {r^,h/h} = 0.14 ± 0.06, for ?i = 2 bulges 
{rc,b/h) = 0.25±0.17,' for n = 3 bulges {r^,b/h) = 0.35±0.27, 
and for ?i = 4 bulges (ro,b//i) = 0.54 ± 0.54. 

Since n is an indicator (in statistical sense) of morpho- 
logical type, the changes of r^^h/h with the shape of the 
bulge implies that Hubble sequence is not scale- free. The 
same conclusion has been recently made by Gadotti (2009) 
for his sample of nearly 1000 galaxies. Moreover, he was 
able to distinguish in his sample two subsamples of classi- 
cal bulges and pseudobulges. Gadotti (2009) have used the 
Kormendy relation (Kormendy 1977) for identifying pseu- 
dobulges and bulges by their loci in the (/ic.b j3 - '"e.b plane. 
The corresponding loci are quite different for these two types 
of bulges. Gadotti (2009) found that re^h/h for galaxies of 
both subsamples was clearly correlated with the bulge-to- 
total luminosity ratio, although the corresponding relations 
for pseudobulges and classical bulges were offset. As the 
bulge-to-total luminosity ratio is related to the morphologi- 
cal (Hubble) type of a galaxy, the existence of such relations 
rule out the hypothesis of a scale-free Hubble sequence. It 
means that the widely discussed secular evolution scenario 
of the bulge formation should be improved to include the 
variation of r^^h/h with n (or other indicators of the type) 
even in the case of pseudobulges. 

The disc flattening can be measured directly only in 
edge-on galaxies. It is commonly expressed as the ratio of 
disc scaleheight zq to disc scalelength h. This ratio shows a 
weak trend with Hubble type (de Grijs 1998) , but there is a 



* The mean effective surface brightness within the circle of effec- 
tive radius b- 



Table 7. The best fitting results for correlations between 
bulge and disc scale parameters 
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fairly clear correlation between both scale parameters (e.g. 
Kregel et al. 2002). The result of Kregel et al. (2002) can be 
transformed into the mean value of (h/zo) = 3.7±1.1 for the 
volume corrected distribution. Our data also demonstrates 
such a correlation for galaxies of all types (Fig. llOb l . It is 
remarkable that the line of the linear regression is drawn 
through the point close to (0,0). 

We used the data of photometric parameters h and zq 
in BM02 {Ks filter) of their reliable subsample and found 
the following equation of linear regression: zq = (0.210 ± 
0.026) h - (0.017 ± 0.258) with correlation coefficient r = 
0.727, which is very close to our result (Table[7]). 

Thus, the galaxies with large disc scalelengths have on 
average larger disc scaleheights and larger effective radii of 
their bulges. However, one should keep in mind the fact that 
the spread in the h — zq plane is large. This scatter is re- 
flected in a wide and remarkably flat distribution of the ratio 
h/zo presented in Fig.[5}l. We return to this question in Sec- 
tion [53] 

As is seen in Fig. llOb the disc scaleheight tends to in- 
crease with the bulge effective radius. The correlation be- 
tween Zo and rc,b is new and was not described earlier. 
We found the following ratios of mean r^b/zo: ('"o,b/-zo) = 
1.23 ± 0.48 for the J-band, {r^,h/zo) = 1.57 ± 0.62 for the 
//-band, (re.b/^o) = 1.49 ± 0.60 for the /r^-band. 

5.3 Well-known relationships 

Here we investigate the general scaling relations for spiral 
galaxies such as the TuUy-Fisher relation, the Fundamental 
Plane (FP) of discs and the Photometric Plane of bulges 
(PhP). In this subsection we study the TuUy-Fisher rela- 
tion and the Fundamental Plane of discs and in the next 
subsection we specially investigate the Photometric Plane 
of bulges. 

5.3.1 The Tully-Fisher relation 

The Tully-Fisher (TF) relation links the luminosity and the 
maximum rotation velocity Urot of galaxies. The small scat- 
ter in the Tully-Fisher relation allows us to use it as a 
method for measurements of distances to spiral galaxies. To 
determine the slope and the zeropoint of the TF relation we 
used the total luminosities of the sample galaxies taken from 
2MASS (Joxt, J/ext or /^oxt). They were derived from the 
isophotal magnitudes by using the iC's-band 20 mag/ arcsec^ 
elliptical radius (J20, H20, K20) and then were extrapolated 
to roughly four times the disc scalelength h (for the detailed 
description see Jarrett et al. 2000 and Cutri et al. 2006). 
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Figure 9. (a) The disc scalclcngth h vs. the bulge effective radius b, (b) ratio vs. Sersic index n in the A's-band. Open squares 

arc related to the whole sample and the filled squares are related to the complete sample. The solid line corresponds to the regression 
line for the complete sample. The errorbars indicate the median error for all galaxies. 




Figure 10. (a) The disc scaleheight zq vs. the disc scale length h, (b) the disc scale height zq vs. the bulge effective radius rg ^ (Ks). 
Symbols as in Fig. [9] 



Unfortunately, measured rotation velocities are not 
known for all the sample galaxies. In the J-band the num- 
ber of galaxies with known rotational velocities is 107, in the 
//-band there are 110 galaxies and in the Ks-hand there are 
116 objects with known velocities. The maximum rotational 
velocity corrected for inclination v^ot was taken from LEDA. 

The absolute magnitude in different bands was calcu- 
lated as 

Ma = rriA - - 25 - 5 Ig(DL) , (10) 

where mx is the total apparent magnitude in each passband 
(Jext, //ext or ifcxt), /^L IS the luminosity distance and A\ 
is the Galactic extinction taken from Schlegel et al. (1998). 
The luminosity is expressed as 

IgLA =0.4(M(^ -Ma), (11) 

where M^ = 3.64 mag, Mg = 3.32 mag and = 3.28 
mag (Binney & Merrifield 1998). We fitted the TF relation 
coefHcients by taking the bisector between the lines of linear 
regression obtained through the standard LSQ and inverse 
LSQ methods. The best fitting results are presented in the 
Table U The Ti's-band TF relation is plotted in Fig.lTTK. 

There are many papers with the TF relation in the NIR 
bands for early- and late-type galaxies. Here we compare our 
results with those of some works based on 2MASS observa- 



tions. The results by Karachentsev et al. (2002), Courteau 
et al. (2007), de Rijcke et al. (2007) and Masters, Springob 
& Huchra (2008) are presented in Table [8] 

As was shown by Masters et al. (2008) the TF relation 
coefficients depend on galaxy morphology. The relationship 
is steeper for later type spirals than earlier type spirals. The 
TF zeropoint is dimmer for late-type galaxies. We did not 
divide our sample on morphological types to get the TF 
relation for various galactic types. Our sample contains an 
almost equal number of early- and late-type galaxies. That 
is why the coefficients found are somewhat different from 
those given in other papers. 

5.3.2 The Fundamental Plane of discs 

Following Moriondo et al. (1999) we fit a relation for disc 
parameters as: 

Ig /i = a Ig Wrot -f fe S'o.d -I- c , (12) 

involving the disc scalelength, the maximum rotational ve- 
locity, and the deprojected central surface brightness So, 
analogous to the FP of elliptical galaxies. To define coef- 
ficients of the linear approximation we used the standard 
least squares method. From the complete subsample we ex- 
cluded 7 galaxies: NGC 7232 and NGC 5775 are interact- 
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Table 8. Comparison between some published 2MASS TF relations IgL = algtv, 



Reference 


J- band 


//-band 


iC-band 


Comments 




a 


b 


a 


b 


a 


b 




This work 


3.43 


2.67 


3.44 


2.83 


3.55 


2.66 


all types 


Karachentsev et al. (2002) 


3.28 




3.48 




3.61 




Based on Jcxt, //cxt 


Courteau et al. (2007) 


2.5 


4.0 


2.6 


4.0 


2.6 


3.9 


64% Sa galaxies 


de Rijcke et al. (2007) 










3.46 


2.44 


early-type galaxies 


Masters et al. (2008) 


3.63 


1.05 


3.61 


1.30 


4.01 


0.36 


Sc galaxies 


Masters et al. (2008) 


3.12 


2.32 


3.12 


2.53 


3.45 


1.76 


Sb galaxies 


Masters et al. (2008) 


2.44 


4.11 


2.43 


4.31 


2.78 


3.48 


Sa galaxies 



and Ks 



RFGC galaxies 



ing, NGC 7814 is a pure SO galaxy with very faint disc, 
NGC 4222 and NGC 100 are low surface brightness galax- 
ies, NGC 5023 and NGC 4183 have very large errors of disc 
parameters. The results of the fitting are given in Table |9l 
The scatter of points is rather large which may be caused by 
observational uncertainties and decomposition errors. The 
coefficients derived in this paper are comparable with those 
reported by Karachentsev (1989) for the /-band: a — 1.4, 
b = 0.28; and by Moriondo et al. (1999) for the //-band: 
a = 1.31±0.19, 6 = 0.25±0.04. Our coefficients as well as the 
coefficients derived by Karachentsev (1989) and Moriondo et 
al. (1999) are not consistent with what would be expected on 
the basis of the virial theorem and a universal mass-to-light 
ratio (a = 2 and b = 0.4). 

We divided galaxies into galaxies with n > 2.2 and 
gala^xies with n < 2.2. In spite of the fact that the shape 
parameter n is related to the bulge model, the fundamental 
planes of discs, which tie only disc parameters and the maxi- 
mum rotational velocity of gas, are different for galaxies with 
different bulges. Perhaps it is due to observational errors and 
errors of decomposition for galaxies of different morpholog- 
ical types (for early-type galaxies rotational velocities and 
parameters of faint discs are less well defined than for late- 
type galaxies). Another reason may be the real difference 
of discs in galaxies with low and high density bulges. This 
effect requires a more careful examination. 

5.4 Photometric Plane and Sersic index 

The Sersic index n represents the parameter quantifying the 
radial concentration of the stellar distribution of bulges. The 
detailed structure of bulges, their formation and evolution- 
ary status is now widely discussed. Here we present a cor- 
relation between bulge absolute magnitude A/b and n. The 
plot is shown in Fig. 112b .. The presented correlation is well- 
known (e.g. Andredakis et al. 1995; Graham 2001; Hunt et 
al. 2004; MoUenhoff & Heidt 2001; MoUenhoflF 2004) and 
seen in all bands — from optical to near-infrared. It shows 
that exponential bulges are fainter than the de Vaucouleurs' 
bulges, although the scatter in the plot is large. 

There are a lot of trends of galactic structural parame- 
ters with the Sersic index of a bulge n, but the most promi- 
nent difference in parameters arises when one divides the 
sample under investigation into the subsample of classical 
bulges and pseudobulges. As was shown above, perhaps, the 
difference between classical bulges and pseudobulges is re- 
flected in the intrinsic flattening and the 3D shape of bulges, 
in the ratio r^^h/h and in building fundamental planes for 
discs. 

The division of bulges into classical bulges and pseu- 



dobulges was described in many works (see Kormendy & 
Kennicutt 2004 for review). Classical bulges appear to be 
similar to elliptical gala^cies and have similar properties. It 
is assumed that these systems were built via minor and ma- 
jor merging. The galaxies with pseudobulges are classified 
by the presence of some features within a bulge, e.g. nuclear 
bars, nuclear spirals and nuclear rings (Erwin & Spark 2002; 
CaroUo et al. 1997; Fisher & Drory 2008). Pseudobulges are 
thought to be formed via disc instabilities and secular evolu- 
tion and have disc-like apparent flattening. One of the main 
differences between classical bulges and pseudobulges is in 
the shape parameter n. Pseudobulges as usual have n < 2 
and classical bulges have n ^ 2 with little-to-no overlap. 

The images of galaxies under investigation have low res- 
olution and inner features inside a bulge are smoothed. That 
is why we were searching for correlations between obtained 
structural parameters of the bulges with different Sersic in- 
dex n, not with morphological features. Here we describe the 
difference between bulges and pseudobulges while building 
the photometric plane for them. 

Khosroshahi et al. (2000a, b) found a tight correlation 
of the Sersic indices n with the central surface brightness 
/^o,b and the effective radius of a bulge r^b- They called this 
relationship the Photometric Plane. 

Following Khosroshahi et al. (2000a, b) we performed 
the least-squared fit which expresses Ig n as a linear combina- 
tion of two other parameters: Ig n = a Ig rc + & /lo.b + c. From 
the complete subsample we exclude galaxies with very small 
bulges and 2 bright fiat galaxies of late types: NGC 4244 
and NGC 4517. 

For the complete sample the best fitting results are pre- 
sented in Table 1101 Khosroshahi et al. (2000a, b) found the 
Photometric Plane for 26 bulges of predominantly early- 
type galaxies in the K-hanA: \gn = (0.130 ± 0.040) Ig ^cb — 
(0.073±0.011)^io,b+(l-21±0.11). MoUenhoff & Heidt (2001) 
obtained for their sample of 41 predominantly late-type spi- 
rals in the JHK bands: Ign = 0.1871grc,b — 0. 081/^0, b + 
1.34 ±0.10 with the correlation coefficient r = 0.91. Finally, 
Mendez-Abreu et al. (2008) gave Ign = 0.17(±0.02) Igre.b - 
0.088(±0.004)^io,b + 1.48(±0.05) for their sample of 148 un- 
barred SO-Sb galaxies in the J-band. Within error bars, the 
differences between our coefficients for all our complete sam- 
ple and those of Khosroshahi et al. (2000a, b) , MoUenhoff & 
Heidt (2001) and Mendez-Abreu et al. (2008) are not signif- 
icant, at least in the K-hanA. 

The surprising result was the non-linear correlation of 
the PhP for galaxies with various values of n (Fig. 112b ). 
From our main sample we extracted the subsample of galax- 
ies with Ign > 0.2 (n > 1.58) and obtained the best-fit Pho- 
tometric Plane (see Table llOp . It occurred that the galax- 



© 2009 RAS, MNRAS 000.[TlfT9l 



14 A. V. Mosenkov, N. Ya. Sotnikova and V. P. Reshetnikov 



Table 9. The TuUy-Fishcr relations and disc fundamental planes for com- 
plete sample galaxies in the J, H, and /^s-bands. The fundamental planes 
are presented for all complete sample galaxies, for galaxies with n > 2.2 
and for galaxies with n < 2.2. 



J: 
H: 
Ks-. 


IgL = (3.43 ± 0.29) Ig Vrot + (2.67 ± 1.39), r = 0.82 
IgL = (3.44 ± 0.28) IgUrot + (2.83 ± 1.32), r = 0.82 
IgL = (3.55 ±0.30) Igiirot + (2.66 ± 1.41), r = 0.81 




For all galaxies of the complete sample with known velocities: 


J: 
H: 
Ks-. 


Igh = (1.017 ± 0.101) IgDrot + (0.173 ± 0.021)So,d - 
Ig/i = (1.044 ± 0.092) \gVrot + (0.190 ± 0.020)S'o,d - 
\gh = (1.077 ± 0.104) Igvrot + (0.202 ± 0.023)5o',d - 


(5.04 ±0.11) 
(5.28 ±0.11) 
(5.53 ±0.12) 


For 


galaxies of the complete sample with n > 2.2: 




J: 
H: 
Ks-. 


Igh = (0.729 ± 0.265) Igtirot + (0.204 ± 0.054)5o,d - 
lg/i= (0.790 ±0.149) \gVrot + (0.193 ± 0.029)So,d - 
\gh = (0.750 ± 0.149) Igvrot + (0.209 ± 0.031)5o'd - 


(4.92 ±0.14) 
(4.70 ±0.11) 
(4.83 ±0.12) 


For 


galaxies of the complete sample with n < 2.2: 




J: 
H: 
Ks-. 


lgh = (1.218 ± 0.122) IgDrot + (0.213 ± 0.027)So,d - 
Ig/i = (1.224 ± 0.149) \gVrot + (0.230 ± 0.025)5o,d - 
\gh = (1.281 ± 0.146) Igvrot + (0.231 ± 0.026)5o',d - 


(6.29 ±0.09) 
(6.45 ± 0.08) 
(6.53 ± 0.08) 




1.8 2 2.2 2.4 2.6 0.2 0.4 0.6 0.8 1 

Ig V„t (km s-1) 0.751g V„t + 0.21So,d - 4.83 



Figure 11. (a) The TuUy-Fisher relation in the i^^-band. Symbols as in Fig. [9] (b) The disc fundamental plane for galaxies with n > 2.2 
in the Ks-hand (filled squares). It is shown that galaxies with n < 2.2 compose a different fundamental plane. These galaxies are plotted 
with open triangles. Only galaxies of the complete sample are shown. 



ies with Ig n > 0.2 lying on the PhP, that is built for these 
galaxies, has a small scatter, but the galaxies with Ig ti < 0.1 
(n < 1.26) do not correspond to this plane and form their 
own plane (see Table [10)) . Inside the window 0.1 < Ig m < 0.2 
there are points which do not correspond to these two planes. 
The regression of the PhP is non-linear, but we approximate 
this dependence as a superposition of two linear regressions. 
The definite difference between the photometric planes of 
the bulges with different shape parameters n is shown in 
Fig. 112b . The error bars are plotted as an indicator of the re- 
ally different location of bulges and pseudobulges in the pho- 
tometric planes. At a first glance the errors in the leftmost 
region of the plot are large, but bulges with very small n 
are clearly recognized in surface brightness profiles of galax- 
ies. What more, the relationship for bulges in the rightmost 
region of the plot begins to deviate from a linear trend for 
values of n significantly above 1. The main difference be- 
tween our sample and the samples of the above-mentioned 
authors is the range of the shape parameters n. Our sam- 
ple contains a substantial amount of bulges with Ign < 0.1. 
For other samples the low boundary for n lies at Ign > 0.2, 
that is why the curvature of the Photometric Plane towards 
small values of n for bulges was not noticed earlier. 



This result is quite new and surprising. We realise that 
the points on the lower-left of the figure have bulges with 
n < 1, which are somewhat controversial and need further 
investigation, which will be done in a future paper. But 
a similar planar relation associated with the Photometric 
Plane was found for bright ellipticals, dwarf ellipticals and 
lenticular galaxies (Khosroshahi et al. 2000a; Khosroshahi 
et al. 2004; Ravikumar et al. 2006; Barway et al. 2009). If 
the range of the Sersic parameter n for ellipticals stretches 
to very small values of n, the curvature is clearly observed 
(Khosroshahi et al. 2004; Ravikumar et al. 2006). 

The origin of the tight relation betveen n and the linear 
combination of Ig and /io,b as well as the curvature of this 
relation are not completely clear. Aceves, Velasquez & Cruz 
(2006) argued that the curvature was related to an intrinsic 
property of a Sersic profile. They wrote an expression for 
the total luminous matter associated with a Sersic profile 
in log-space, involving n, Igro and /io,b, and showed that 
for a set of galaxies with equal luminosities there were non- 
constant terms in this expression. These terms introduce a 
systematic change in a PhP-like expression. 

It is not clear whether two distinct parts of a curve in 
Fig. 112b reflect the different origin of bulges lying in the cor- 
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responding regions. However, the coUisionless merger rem- 
nants of disc galaxies obtained in A^-body simulations and 
fitted by a Sersic profile fairly well reproduce the slope of 
the PhP in the region of Ig n > 0.2 (Aceves et al. 2006). We 
put on the data of Aceves et al. (2006) in our Fig. 112b and 
found the numerical data are consistent with our observa- 
tional ones. 

The curvature of the PhP towards small values of n may 
reflect the quite different nature of such bulges, formed, for 
example, via secular evolution of discs. One of the signs of 
secular evolution is the presence of a bar or other structures 
(such as rings, double exponential discs). We have made the 
assumptions about the existence of a bar or other struc- 
tures in our sample galaxies on the basis of the ID surface 
brightness profiles along the major axis and on the residual 
images where additional components, that were not included 
in the decomposition, become apparent. The galaxies with 
or without bars are represented among the galaxies with 
Ign > 0.2 (classical bulges) as well as among the galaxies 
with Ig n < 0.2 (pseudobulges). So the shape of the PhP may 
describe the deeper nature of these objects. Marquez et al. 
(2001) have suggested a theoretical explanation of the PhP 
of ellipticals. They calculated the specific entropy of ellipti- 
cal galaxies (Entropic Surface), took into account a scaling 
relation between the potential energy and mass (Energy- 
Mass Surface) and showed that the PhP arose as an inter- 
section line of these two planes demonstrating a curvature 
consistent with the observational data in the limited range 
of values that has been considered for n. All these arguments 
may be valid for bulges of all types too. 

5.5 The dark halo and disc flattening 

In spite of a good correlation between the disc scalelength h 
and the disc scaleheight zq, the scatter of the ratio h/zo is 
fairly large. The ratio h/ zo characterizes a relative thickness 
of a galaxy, or its flattening. Fig. [5jl shows that h/zo varies 
from ~ 8 for thin galaxies to 1 — 2 for plump discs. Disc 
flattening appears to correlate with the morphological type 
(e.g. de Grijs 1998) and the global HI content (Reshetnikov 
& Combes 1997). Late-type and gas-rich galaxies are, on av- 
erage, thinner than early- type systems, but this correlation 
is rather weak. Although the discs seem to be larger and 
thicker with a higher maximum rotation (van der Kruit & 
de Grijs 1999; Kregel et al. 2002, 2005), the ratio h/zo is 
not correlated with the rotation velocity (Zasov et al. 2002; 
Kregel et al. 2002, 2005). 

Zasov et al. (2002) and Kregel et al. (2005) found that 
the disc flattening shows a definite trend with the ratio of 
dynamical mass to disc luminosity provided the dynami- 
cal mass is defined as the total mass enclosed within the 
sphere of radius that is equal to four disc scalelengths: 
Mtot = 4/iWrot/G. The total mass includes the mass of the 
disc Md, the dark halo Mh and the bulge Mb. So, the ratio 
of dynamical mass to disc luminosity determines the relative 
mass of a spherical component or the relative mass of the 
dark halo for bulgeless galaxies. The found correlation im- 
plies that relatively thinner discs in bulgeless galaxies tend 
to be embedded in more massive dark haloes. 

Kregel et al. (2005) argued that the analytical collapse 
model of disc galaxy formation (White & Rees 1978; Fall 
& Efstathiou 1980; Blumenthal 1986; Mo, Mao & White 



1998), the local Toomre's (1964) stability criterion for a disc 
and the suggestion about vertical dynamical equilibrium, 
provided the constant value of the ratio of the vertical to 
the radial velocity dispersion (Jz/ctr, lead to the correlation 
between the disc flattening and the ratio of dynamical mass 
to disc luminosity. 

Zasov et al. (2002) suggested that discs are marginally 
stable against the growth of perturbations in their planes 
and bending perturbations. This stability against bending 
perturbations means that the ratio CTz/ctr must be larger 
than some analytical threshold: a^ja-R. > 0.3 — 0.4 (Poly- 
achenko & Shukhman 1977; Araki 1985). A''-body simula- 
tions give a somewhat larger value for this threshold — 
o-z/tJR ~ 0.6 - 0.8 (e.g. Sotnikova & Rodionov 2003, 2006). 
For a disc in vertical dynamical equilibrium Zasov, Makarov 
& Mikhailova (1991) and Zasov et al. (2002) deduced a sim- 
ple expression for hj zq. We rewrite it in the following form: 

h 1 Mtot /a Mtot 

^0 \Q{G^lo^)f Md [Q(a./aR)]2 id ' 

where Q is Toomre's (1964) stability parameter, and /a is 
the mass-to-luminosity ratio in a flxed band. The linear rela- 
tion between the ratio /i/zq and the ratio of dynamical mass 
(or dark halo mass) to disc luminosity arises for marginally 
stable discs with Q and a-z./on near their thresholds. 

To verify a linear trend (|13p . obtained analytically, 
Zasov et al. (2002) studied two different samples of edge- 
on bulgeless galaxies with known structural parameters in 
the R and Ks bands. They used the HI line width W50 as 
a measure of dynamical mass and concluded that discs be- 
come thinner with increasing mass fraction of their dark 
halos. For our sample, which contains a great many galaxies 
with massive bulges, we performed the same analysis. 

Fig. 113b , shows the distribution of the sample galaxies 
over the ratio of the dynamical mass to the stellar mass 
M, = Md -I- Mb. We adopted mass-to-light ratios fj = 1.5, 
/h = 1.0, and /x, = 0.8 Mq/Lq (McGaugh et al. 2000). 
These mass-to-light ratios are consistent with the maximum 
disk best fitting the rotation curves for bright galaxies (e.g. 
Bottema 1999). There is a hint of bimodality in Fig |13b ,. This 
may reflect the presence of two different families of galaxies 
with different bulges in our sample. 

The distribution of the ratio Mtot/M, ranges over the 
values from ~ 1 to ~ 8 — 10. The distribution of the ratio 
h/zo is varying in the same range (Fig. [5}1). As the rela- 
tion (|13l) deals not only with the contribution of the dark 
halo to the total mass but with the relative mass of spherical 
components, including a bulge, we choose the ratio Mtot /Md 
to demonstrate a trend for the disc flattening with the rela- 
tive mass of a spherical component fFig. 113b and Table [TT|l . 
We excluded 6 points from fltting of Mtot/Md vs. h/zO: 
NGC 7814, NGC 4222, NGC 5965, NGC 4183 and NGC 100, 
that have faint disc, and NGC 5981 with a prominent bar. 

The scatter of points in Fig. 113b is rather large. It may 
be caused by uncertainties in velocities and mass-to-light ra- 
tio determinations, decomposition errors and disc relaxation 
processes, that can thicken a disc just above the threshold 
for marginal stability. But some of this scatter must be real. 
Sotnikova & Rodionov (2005) concluded that the presence 
of a compact bulge is enough to suppress the bending insta- 
bility that leads to the disc thickness increasing. A series of 
A'^-body simulations with the same total mass of a spherical 
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Table 10. The Photometric Planes for bulges of the complete sample 



For all galaxies of complete sample 

J: \gn= (0.053 ±0.059) lgre,b - (0.101 ± 0.004)/io,b + (1.63 ±0.06) 
H: Ign = (0.106 ± 0.047) Igr^.b - (0.084 ± 0.004)Ato,b + (1.37 ±0.05) 
Ks-. Ign = (0.160 ± 0.058) lgrc',b - (0.084 ± 0.005)Ato'b + (1.34 ±0.07) 



Ign ^ 0.2 

J: Ign = (0.030 ±0.018) lgre,b 
H: Ign = (0.018 ±0.014) lgrc,b - 
Ks:lgn = (0.022 ± 0.013) Ig r^'^b 



' (0.067 ± 0.002)/io,b + (1-23 ± 0.03) 
(0.062 ± 0.002)Ato,b + (1-14 ± 0.02) 
' (0.061 ± 0.002)Ato',b + (1.11 ± 0.02) 



Ign < 0.1 

J: Ign = (0.137 ±0.134) lgre,b 
H: Ign = (0.201 ± 0.162) lgre,b ■ 
Ks:lgn = (0.467 ± 0.214) Igro'b 



- (0.142 ±0.039)/io,b - 
(0.172 ± 0.051)Ato,b - 
■ (0.204 ± 0.075)Ato'b ■ 



' (2.31 ±0.67) 
(2.76 ± 0.83) 
■ (3.21 ± 1.19) 
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Figure 12. (a) Bulge Scrsic index n vs. bulge absolute magnitude in the i^s-band. Symbols as in Fig. [9] (b) The Photometric Plane 
built for bulges with Ign ^ 0.2 of the complete subsample in the Ks-band. The points corresponding to bulges with Ign < 0.1 do not 
lie in this plane. The model merger remnants of disc galaxies are shown with open circles and triangles (Aceves at al. 2006). The model 
data are shifted by a constant along the horizontal axis. 



Table 11. Correlations between zq/H and Mtot/Mj 



J: 


h/zo 


= (0.634 ± 0.100)Mtot/Md - 


h (1.63 ± 1.16), r 


= 0.58 


H: 


h/zo 


= (0.696 ± 0.120)Mtot/Md H 


- (1.76 ± 1.51), r 


= 0.51 


Ks 


h/zo 


= (0.756 ± 0.135)Mtot/Md - 


\- (1.79 ± 1.92), r 


= 0.46 



component (dark halo + bulge) were performed. The final 
disc thickness was found to be much smaller in the simula- 
tions, where a dense bulge is present, than in the simulations 
with bulgeless systems. The results of A^'-body simulations 
of discs starting from an unstable state were summarized by 
Sotnikova & Rodionov (2006). They plotted the ratio zo/h 
versus (Mb + Mb)/Md and showed that there is a clear scat- 
ter in this relation, in spite of the same model mass for a 
spherical component (Mb + Mb). We are planning to com- 
pare in details our observational data with the results of 
simulations in a forthcoming paper. 



6 SUMMARY AND CONCLUSIONS 

We constructed the largest sample of edge-on spiral galax- 
ies with performed 2D bulge/disc decomposition of 2MASS 



galaxies images in J-, H- and Jfa-passband using the 
BUDDA v2.1 package. We extracted global bulge and disc 
parameters for all objects. This first paper describes our 
edge-on galaxy sample. The main point of this work is in 
the first carefully fitting of edge-on spiral galaxies of early 
and late morphological types in 3 near-infrared bands. 

Edge-on galaxies are of great interest because they pro- 
vide a unique possibility to obtain information about the 
vertical structure of a disk and a bulge. The disc scaleheight 
zo together with the disc scalelength h determine the relative 
thickness of a stellar disc for each galaxy. This ratio as well 
as the bulge flattening may put constraints on the disc 
and bulge formation processes and their secular evolution. 

The results of the sample parameter analysis, that are 
specific for edge-on galaxies, can be summarized as follows: 

(i) Our data demonstrate a clear correlation between the 
disc scaleheight zo and the disc scalelength h for all type 
galaxies. The distributions of the sample galaxies over the 
ratio of h/zo have a median value of 3.5 for J-band and 
3.9 for H and Ks-haxids. But there is a substantial scat- 
ter of this ratio (the distribution is fairly flat). The ratio 
h/zo stretches over a large range — from very thin galaxies 
to plump ones. There is a correlation between the relative 
thickness of stellar discs h/zo and the relative mass of a 
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spherical component, including a dark halo. This correla- 
tion was known previously for bulgeless galaxies (Zasov et 
al. 2002) or samples with predominantly bulgeless galaxies 
(Kregel et al. 2005) and was argued to arise from marginally 
stable discs. Our sample is much larger than the samples of 
Zasov et al. (2002) and Kregel et al. (2005) and more reli- 
able for statistical analysis. Here, we confirm the correlation 
under discussion. What is more, we do this not for bulgeless 
galaxies but for galaxies with massive bulges. We conclude 
that the total mass contained in spherical components may 
be one of the factors that determines the final steady state 
disc thickness. 

(ii) The bulge flattening in the vertical direction gb di- 
vides the sample into two different families — triaxial, nearly 
prolate bulges and close to oblate bulges with moderate flat- 
tening. The Sersic index threshold n ~ 2 can be used to 
identify these two bulge types. 

(iii) We found a correlation between zo and ro,b. It is new 
and has not been described previously. It means that the 
disc flattening is linked with the bulge structure. 

Many of our results are in good agreement with the 
results of other authors, but we found several new relations. 

(i) In accordance with previous studies we found that the 
scale parameters of bulges and discs are tightly correlated. 
The disc scalelength linearly increases when the bulge ef- 
fective radius increases. However, there is a clear trend for 
the ratio r^^i^lh to increase with n. As n is an indicator of 
the Hubble type, such a trend unambiguously rules out the 
widely discussed hypothesis of a scale-free Hubble sequence. 

(ii) There is a hint that the fundamental planes of discs, 
which link only disc parameters and the maximum rota- 
tional velocity of gas, are different for galaxies with differ- 
ent bulges. This may indicate a real difference of discs in 
galaxies with low and high density bulges. 

(iii) The most surprising result arises from the investiga- 
tion of the Photometric Plane of sample bulges. The bulges 
with 71 > 2 populate a narrow strip in their Photometric 
Plane. However, there is a difference in behavior of this plane 
for bulges with n > 2 and n < 2. The plane is not flat and 
has a prominent curvature towards small values of n. For 
bulges this fact was not noticed earlier. This result may be 



due to the physical distinction between classical bulges and 
pseudobulges. As the Photometric Plane is proposed to be 
used as a distance indicator, the correct shape of the plane 
is of no little interest. 

Some of results were described only briefly, but it is clear 
that our sample is very useful for further detailed studying 
and modelling of edge-on spiral galaxies. In the forthcom- 
ing papers we are planning to discuss properties of edge-on 
galaxies in more detail in the context of their bulge and disc 
formation and joint evolution. 
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Table 1. The Sample. 



# 


Galaxy 


2MFGC 


Type 


D 

(Mpc; 




(1) 


(2) 


(3) 


(4) 


1 


UGC 529 


633 


SO-a 


70.4 


2 


NGC 504 


1059 


SO 


53.2 


3 


NGC 527 


1070 


SO-a 


74.7 


4 


NGC 565 


1133 


Sa 


57.0 


5 


UGC 1120 


1193 


Sab 


58.8 


6 


ESQ 353-G020 


1202 


SO-a 


61.7 


7 


UGC 1166 


1243 


SO-a 


60.8 


8 


ESQ 013-G024 


1422 


Sab 


— 


9 


UGC 1531 


1557 


Sb 


105.0 


10 


IC 207 


1660 


SO-a 


62.1 


11 


NGC 861 


1758 


Sb 


106.0 


12 


UGC 1938 


1930 


Sbc 


82.0 


13 


NGC 955 


1965 


Sab 


17.0 


14 


NGC 960 


1984 


Sb 


63.1 


15 


MCG-02-07-038 


2044 


SO 


62.2 


16 


NGC 1029 


2102 


SO-a 


46.0 


17 


UGC 2304 


2260 


Sbc 


93.7 


18 


MCG-02-09-013 


2605 


SO-a 


61.2 


19 


ESQ 005-G004 


4913 


Sb 


26.1 


20 


NGC 1381 


2983 


SO 


22.2 


21 


NGC 1394 


3025 


SO 


55.5 


22 


NGC 1401 


3030 


SO 


18.8 


23 


NGC 4256 


9696 


Sb 


35.7 


24 


MCG-02-08-012 


2274 


Sc 


68.8 


25 


NGC 1529 


3366 


SO 


68.7 


26 


UGC 7388 


— 


SBb 


89.8 


27 


ESQ 251-G028 


3869 


Sb 


158.0 


28 


NGC 4866 


10288 


SO-a 


31.2 


29 


ESQ 555-G023 


4860 


SO-a 


— 


30 


ESQ 555-G032 


4905 


Sb 


— 


31 


ESQ 087-G004 


5097 


Sbc 


118.0 


32 


UGC 4012 


6163 


Sb 


130.0 


33 


UGC 4190 


6390 


SO 


69.9 


34 


UGC 4507 


6803 


Sb 


141.0 


35 


UGC 4526 


6822 


Sab 


62.4 


36 


ESQ 433-G012 


7223 


Sbc 


115.0 


37 


ESQ 373-G013 


7473 


SO-a 


40.7 


38 


NGC 2946 


7484 


SBb 


122.0 


39 


UGC 5345 


7714 


Sb 


53.5 


40 


IC 615 


8108 


Sb 


133.0 


41 


ESQ 501-G056 


8260 


SO 


53.0 


42 


NGC 3431 


8449 


SABb 


75.3 


43 


MCG-01-28-005 


8496 


Sab 


91.0 


44 


CGCG 267-039 


8658 


Sb 


88.2 


45 


NGC 3539 


8691 


SO-a 


132.0 


46 


NGC 1184 


2685 


SO-a 


30.8 


47 


NGC 4149 


9578 


Sb 


43.4 


48 


MCG-02-32-006 


9725 


SBc 


60.0 


49 


NGC 4686 


10067 


Sa 


69.5 


50 


MCG-02-33-076 


10273 


Sab 


47.0 


51 


UGC 8119 


10315 


Sbc 


117.0 


52 


IC 844 


10368 


SO 


43.6 


53 


MCG-01-33-076 


10374 


Sc 


58.4 


54 


IC 5096 


16137 


Sbc 


41.2 


55 


IC 881 


10662 


Sa 


96.9 


56 


NGC 5119 


10737 


SO-a 


43.6 


57 


UGC 8462 


10813 


Sab 


132.0 


58 


UGC 8704 


11102 


Sa 


140.0 


59 


NGC 5290 


11096 


Sbc 


37.4 


60 


NGC 5553 


11621 


Sa 


63.8 
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Table 1. continued. 



# 


Galaxy 


2MFGC 


Type 


D 

(Mpc) 




(1) 


(2) 


(3) 


(4) 


61 


ESO 384-G014 


11249 


SO 


57.2 


62 


IC 1110 


12300 


Sa 


46.0 


63 


ESO 387-G029 


12429 


Sab 


60.6 


64 


UGC 11060 


14135 


Sa 


61.1 


65 


NGC 4235 


9679 


Sa 


37.4 


66 


UGC 11614 


15682 


Sb 


105.0 


67 


IC 5181 


16757 


SO 


24.1 


68 


NGC 7232 


16791 


SO-a 


23.1 


69 


ESO 471-G024 


17877 


SO-a 


113.0 


70 


NGC 482 


1009 


Sab 


85.0 


71 


NGC 522 


1102 


Sbc 


32.8 


72 


NGC 585 


1169 


Sa 


68.9 


73 


NGC 653 


1286 


Sab 


69.9 


74 


NGC 684 


1382 


Sb 


44.1 


75 


NGC 1596 


3625 


SO 


20.5 


76 


IC 2085 


3680 


SO-a 


13.4 


77 


NGC 2295 


5413 


Sab 


26.9 


78 


ESO 311-G012 


6172 


SO-a 


18.2 


79 


NGC 3203 


7996 


SO-a 


37.5 


80 


NGC 3404 


8433 


SBab 


66.8 


81 


NGC 5047 


10588 


SO 


88.7 


82 


ESO 321-GOlO 


9597 


Sa 


46.5 


83 


ESO 416-G025 


2246 


Sb 


65.0 


84 


NGC 3390 


8406 


Sb 


45.7 


85 


IC 127 


1149 


Sb 


23.1 


86 


NGC 3628 


— 


Sb 


16.2 


87 


NGC 3717 


— 


Sb 


28.2 


88 


NGC 5965 


12550 


Sb 


46.8 


89 


NGC 4013 


9412 


Sb 


14.5 


90 


ESO 446-G018 


11466 


Sb 


67.5 


91 


ESO 512-G012 


11898 


Sb 


49.5 


92 


NGC 4710 


10109 


SO-a 


19.6 


93 


NGC 7814 


— 


Sab 


9.5 


94 


ESO 240-GOll 




Sc 


35.7 


95 


NGC 4565 


— 


Sb 


20.7 


96 


NGC 5746 




SABb 


26.4 


97 


NGC 891 




Sb 


4.3 


98 


NGC 5775 


12067 


SBc 


25.7 


99 


NGC 4217 


9661 


Sb 


16.8 


100 


NGC 5908 




Sb 


45.6 


101 


NGC 4244 




Sc 


6.8 


102 


NGC 4570 


9935 


SO 


28.1 


103 


NGC 5907 


12346 


Sc 


10.0 


104 


NGC 3564 


— 


SO 


42.5 


105 


NGC 4266 


9707 


SBa 


25.6 


106 


IC 4202 


10456 


Sbc 


98.5 


107 


NGC 4703 


10103 


Sb 


64.5 


108 


ESO 443-G042 


10371 


Sb 


43.5 


109 


NGC 5308 


11118 


SO 


29.1 


110 


NGC 5365A 


11257 


SBb 


40.5 


111 


IC 1048 


11928 


Sb 


25.2 


112 


UGC 243 


293 


Sb 


66.6 


113 


tibU olo-GOzO 


7524 


bb 


68.9 


114 


NGC 669 


1340 


Sab 


59.5 


115 


NGC 973 


2024 


Sb 


62.4 


116 


IC 335 


2973 


SO 


20.7 


117 


IC 1970 


2982 


Sb 


15.9 


118 


NGC 1380A 


2985 


SO 


20.0 


119 


ESO 362-GOll 


4306 


Sbc 


18.8 


120 


NGC 1886 


4383 


Sbc 


24.3 



Table 1. continued. 


# 


Galaxy 


2MFGC 


Type 


D 










(Mpc) 




(1) 


(2) 


(3) 


(4) 


121 


NGC 2310 


5487 


SO 


17.9 


122 


NGC 2654 


6926 


SBab 


19.8 


123 


NGC 2862 


7315 


SBbc 


58.9 


124 


NGC 3126 


7847 


Sb 


73.2 


125 


NGC 3692 


8974 


Sb 


28.3 


126 


NGC 3957 


9342 


SO-a 


27.1 


127 


NGC 4417 


9805 


SO 


15.8 


128 


NGC 1247 


2619 


Sbc 


50.9 


129 


NGC 5981 


12588 


Sbc 


24.5 


130 


NGC 4026 


9425 


SO 


15.3 


131 


NGC 3501 


— 


Sc 


20.0 


132 


NGC 5023 


10525 


Sc 


8.3 


133 


NGC 360 


761 


Sbc 


30.1 


134 


UGC 1817 


1825 


Scd 


47.2 


135 


MCG-02-10-009 


3079 


Sc 


27.7 


136 


UGC 3326 


4605 


Sc 


55.1 


137 


UGC 3474 


5237 


Sc 


49.1 


138 


NGC 3044 


7660 


SBc 


22.3 


139 


NGC 3279 


8207 


Sc 


23.7 


140 


NGC 4222 


9670 


Sc 


7.7 


141 


ESO 533-G004 


16767 


Sc 


31.2 


142 


NGC 4330 


9747 


Sc 


25.8 


143 


ESO 564-G027 


7159 


Sc 


33.8 


144 


ESO 531-G022 


16371 


Sbc 


42.1 


145 


IC 4393 


11612 


Sc 


40.5 


146 


ESO 263-G015 


7900 


Sc 


38.0 


147 


NGC 4835A 


10246 


Sc 


49.3 


148 


ESO 288-G025 


16587 


Sbc 


30.7 


149 


IC 2058 


3483 


Scd 


18.6 


150 


ESO 201-G022 


3389 


Sc 


55.4 


151 


IC 4871 


15015 


SABc 


25.1 


152 


IC 4484 


11988 


Sc 


61.3 


153 


ESO 340-G008 


15435 


Sc 


36.1 


154 


UGC 12533 


17534 


Sb 


75.8 


155 


ESO 509-G019 


10818 


Sbc 


144.0 


156 


ESO 506-G002 


9712 


Sbc 


57.9 


157 


NGC 4183 


9620 


Sc 


15.8 


158 


NGC 4517 


9881 


Sc 


20.1 


159 


NGC 2357 


5811 


Sc 


33.0 


160 


UGC 5173 


7514 


Sb 


87.6 


161 


NGC 100 


282 


Sc 


6.7 


162 


ESO 121-G006 


4933 


Sc 


17.6 


163 


NGC 2424 


6083 


SBb 


47.3 


164 


IC 2207 


6207 


Sc 


66.7 


165 


ESO 563-G014 


6848 


SBcd 


27.1 


166 


UGC 6012 


8487 


Sbc 


88.8 


167 


ESO 507-G007 


10056 


Sbc 


75.2 


168 


IC 3799 


10098 


Scd 


54.4 


169 


MCG-03-34-041 


10613 


Sc 


40.2 


170 


NGC 5073 


10645 


SBc 


41.4 


171 


IC 4351 




Sb 


39.8 


172 


NGC 5529 


11577 


Sc 


41.4 


173 


NGC 5714 


11872 


Sc 


32.1 


174 


ESO 340-G026 


15509 


Sc 


71.2 


175 


UGCA 150 


7144 


SABb 


29.3 



Columns: (1) Name (the first in NED). (2) 
Number in 2MFGC (2MASS Flat Galaxy Cat- 
alog. (3) Morphological type taken from LEDA 
(Lyon/MeudonExtragalactic Database). (4) Angular- 
size distance [Hq = 73 (km sec~-'^)/Mpc, Cm = 0.27, 
Qa = 0.73). 
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